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ABSTRACT 

The very short and bright flare of 3C 279 detected with Fermi-LAT in 2013 December is tested by 
a model with stochastic electron acceleration by turbulences. Our time-dependent simulation shows 
that the very hard spectrum and asymmetric light curve are successfully reproduced by changing 
only the magnetic field from the value in the steady period. The maximum energy of electrons 
drastically grows with the decrease of the magnetic field, which yields a hard photon spectrum as 
observed. Rapid cooling due to the inverse-Compton scattering with the external photons reproduces 
the decaying feature of the light curve. The inferred energy density of the magnetic field is much less 
than the electron and photon energy densities. The low magnetic field and short variability timescale 
are unfavorable for the jet acceleration model from the gradual Poynting flux dissipation. 

Subject headings: acceleration of particles — quasars: individual (3C 279) — radiation mechanisms: 
non-thermal — turbulence 


1. INTRODUCTION 


Multi-wavelength light curves of blazar flares show 
complex and diversified features. While in some cases 
there is a ti me lag between gamma-r a y and X-ray/optica l 


flares (e.g. 

Blaze 

owski et al.l 12005 

: Fossati et al.l 120081: 

lAbdo et al.l 

2010a 

; iHavashida et al. 

120121). in other cases 


an o rphan flare in a certain wavelength was d etected 
(e.g. iKrawczvnski et al.l[200l: lAbdo et al.ll^lObl) . Even 
if a time-dependent model is adopted, such a vari¬ 
ety of behaviors may be difficult to reproduce with 
by a one-zone model (iKusu no se. Takahya fc L i 1 120001: 
Krawczvnski. Coppi fc Aharonian I 120021: I Asano et all 
20I4D . While spatial gradie nts of the physical parame¬ 
ters in the emission regions (l.laniak et al.ll201^ may ex¬ 
plain some fraction of the lags, some flares have spectral 
evolutions thatare too complex to be modeled, even with 
time-dependent m ulti-zone radiative transfer simulations 
(jChen et al.ll20Ill l. This may imply that inhomogeneous 
emission regions evolve with a longer timescale than a 
dynamical one. Such nontrivial properties in a blazar 
flare make it difficult to probe physical processes such as 
electron acceleration or cooling. 

In 2013 December, the Fermi-Laige Area Telescope 
(LAT) detected one of the most intense flares in 
the gamma-ray band from flat spectrum radio quasar 
(FSRQ) 3C 279, reaching ^ 1 x 1 0~^ ph cm~^ s~^ for the 
integrated flux above 100 MeV (|Havashida et al.l I2015L 
hereafter H15). The flux level is comparable to the his¬ 
torical ma ximum of this source observed at the gamma- 
ray band (jWehrle et al.l 119981 1. The gamma-ray flare 
showed a very rapid variability with an asymmetric time 
profile with a shorter rising time of ~ 2 hr and a longer 
falling time of ^ 7 hr. We can expect that this extraordi¬ 
nary flare was emitted from a sufficiently compact re¬ 
gion that can be regarded as homogeneous, which is 
different from other usual flares. In this case, the de¬ 
caying timescale may directly correspond to the cooling 
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timescale, and the flare is an ideal target for discussing 
the physical processes. 

Another important property of the flare event of 3C 
279, a very hard photon index of T.^ = 1.7 ± 0.1, was 
observed in the > 100 MeV band by Fermi-LAT. Such a 
hard photon index has been rarely observed in FSRQs, 
whose luminosity peak from inverse-Compton (IC) scat¬ 
tering is usually located below 100 MeV. While the mean 
of t he r.v distribution in F SRQs corresponds to about 
2.4 (|Ackermann et ani2015|l . hard photon indices T.^ < 2 
only have been occasionally obs erved in some bright FS¬ 
RQs during rapid flaring events (|Pacciani et akllMll) . In 
order to reproduce the hard photon index by IC scatter¬ 
ing in the fast cooling regime, the index of parent elec¬ 
trons should be much harder than two, which can hardly 
be generated in a normal shock acceleration process. In 
addition, the flare event of 3C 279 indicates a high Comp¬ 
ton dominance parameter L-^/Lsyn > 300, leading to ex¬ 
tremely low jet magnetization with < 10“^ (H15). 

To explain the flare event of 3C 279, rather than 
assuming prompt electron injection by the shock 
acceleration, we propose the stochastic acceleration 
(SA) model, which is phenomenologically equiva¬ 
lent to the second-order Fermi acceleration (Fermi- 
H; e . g. iKatarzvhski et al.l 1200^ iStawarz fc PetrosianI 
I2008I iLefa. Rieger fc Aharonianri20IlL and references 
there in). The SA may be d riven by magnetic reconnec¬ 
tion (jLazarian et al. Il2012[ l. Otherwise, hydrodynami- 
cal turbulences that drive the acceleration are possibly 
induced via t he Kelvin-Helmholtz instability as an ax¬ 
ial mode (e.g. lMizuno. Hardee fc Nishikawall2007D . or the 
Rayleigh-Tay lor and Richtmyer-Meshkov in stabilities as 
radial modes ([Matsumoto fc Masada ll2?)T^ . Broadband 
spectra of blazars in the steady state have been success¬ 
fully fitted with recent SA models (lAsano et all 12014 
iDiltz fc B6ttcherll20l4 iKakuwa et al. II2015D . The flare 
state should be also tested with such models to show the 
wide-range applicability of the SA. This is the first at¬ 
tempt to apply a Fermi-H model to explain both broad- 
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band spectra and light curves of FSRQs simultaneously. 

In this Letter, we perform time-dependent simulations 
of the emissions from 3C 279 with the SA. Starting from 
modeling a steady emission, the gamma-ray flare is re¬ 
produced by decreasing the magnetic held for the steady 
model. We demonstrate that the SA model agrees with 
the observed spectrum and light curve. 

2. NUMERICAL METHODS 

We adopt the numerical simulation code used in 
lAsano et al.l (|2014ll . In this model, a conical outflow with 
an opening angle = l/L, where T is its bulk Lorentz 
factor, is ejected at radius R = Rq from the central en¬ 
gine. The evolutions of electron and photon energy dis¬ 
tributions are calculated in the comoving frame taking 
into account the SA, synchrotron emission, IC scattering 
with the Klein-Nishina effect, adiabatic cooling, 77 pair 
production, synchrotron self-absorption, and photon es¬ 
cape. The SA is charact erized by the energ y diffusion 
coefficient, D{ee) = Kel- lAsano et al.l ^20l4 ) conserva¬ 
tively assumed the Kolmogorov-like diffusion as q = 5/3. 
Here, we adopt, however, q = 2, which corresponds to 
the hard-sphere scattering. This choice leads to reason¬ 
able spectra of 3C 279 without complicated assumptions 
such as nontrivial temporal evolution of the diffusion co¬ 
efficient or electron injection rate. If the cascade of the 
turbulence stops at a certain length scale larger than the 
gyro-radius of the highest-energy electrons, the mean free 
path of electrons becomes comparable to this scale inde¬ 
pendently of electron energies. In this case, the energy 
diffusion can be approximated as the hard-sp here scatter¬ 
ing le.g. iBeresnvak. Yan fc Lazarian Il 2 nilh . The recent 
magnetohydrodynamical simulations accompanying the 
inverse c ascad e shows the k~^ spectrum in turbulences 
(iZrake I 12014 iBrandenburg. Kahniashvili fc Tevzadze I 
I2OI5II . which also support the hard-sphere approxima¬ 
tion. 

The volume we consider is a conical shell with a 
constant width of W = Ro/T, then the isotropically 
equivalent volume V = 47 ri?^IT' (the actual volume 
is Tr9?R^W'). Hereafter, the values in the shell frame 
are denoted with prime characters. During the dynam¬ 
ical timescale W'/c in the plasma frame, electrons are 
injected at a constant rate A/ in the volume V mo- 
noenergetically ( 7 e = 10 ) and acc elerated with a con - 
stant coefficient K'. As done in lAsano et al.l (l2014f) . 
we can consider the temporal evolutions of the injection 
rate and the diffusion coefficient. However, this simple 
model with constant Ad and K' is sufficient to repro¬ 
duce the photon spectrum of 3C 279. The average mag¬ 
netic field in the comoving frame is assumed to behave 
as B' = Bo{R/Ro)~^- The evolution of the photon spec¬ 
trum for observers is computed taking into account the 
relativistic motion and curvature of the jet surface. 

3. STEADY MODEL FOR THE ACTIVE PERIOD IN 2009 

As a reference case, we consider one of the most ac¬ 
tive periods in the gamma-ray band du ring the first two 
years of the Fermz-LAT observations. In iHavashida et id] 
(|2012f) . this period is denoted as period “D” in 2009. Al¬ 
though this period corresponds to an event with a promi¬ 
nent flare, the broadband spectrum in the paper is av¬ 
eraged over five days. If the emission zone is inside the 


broad emission region as suggested by the short variabil¬ 
ity timescale reported in H15, the flare state is signifi¬ 
cantly longer than the variability timescale. Therefore, 
we adopt a steady emission model for period D in 2009. 

By assuming continuous steady ejection of the shells 
from R = Rq, we model the steady photon spectrum, 
though the plasma and its emission evolves with R in the 
shell frame. The model parameters are Rq = 0.023 pc, 
r = 15, AT' = 9 X 10“® s-^ (<acc = l/(2Ar') = O.SSW'/c), 
= 7.8 X 10^® s-i = 0.26{R/Ro)-^ cm-^ s"!), 
and Bq = 7 G. We a dopt the same model as that of 
iHavashida et al.l (|2012f) for the external radiation of the 
broad emission lines with the photon temperature T/jy = 
lOr eV and the energy density U{jy = 8(r/15)® erg cm“^ 
in the shell frame. The cooling timescale in this external 
radiation field is written as 


^cool - 


3meC 

4(JT7eC^UV 



As shown in Figured] (a), electrons are continuously ac¬ 
celerated between R = Rq and 2Rq, then they are rapidly 
cooled via IC scattering after the shutdown of the accel¬ 
eration. The electron spectrum at 2Rq in the low-energy 
region is consistent with the assumed p ower-law index of 
j p = 1 in the broken power-law model of lHavashida et al.l 
(l2012fl . In the highest-energy region, though the Klein- 
Nishina effect suppresses the IC cooling effect, the syn¬ 
chrotron cooling {Uq = 1.9(i?/i?o)“® erg cm“^) prevent 
the acceleration above 100 MeV. The resultant photon 
spectrum well reproduces the observed spectrum from 
far-infrared to gamma-ray bands (see Figure[T](b)). The 
low-energy cutoff at ~ 0.02 eV is due to the synchrotron 
self-absorption. The X-ray flux is originated from the 
synchrotron self-Compton (SSC) emission. Those X-ray 
data strongly constrain the emission radius Rq. 

Thus, our SA model can naturally produce a hard elec¬ 
tron spectrum, and the steady photon spectrum agrees 
with the observed one in 2009. Based on this result, we 
will probe the intensive flare in 2013 in the next section. 


4. FLARE MODEL IN 2013 

The most intensive flare denoted as period “B” (on 
MJD 56646) in HI 5 shows a very hard spectrum with 
r.y ~ 1.7 and a short variability with an hourly scale 
in the gamma-ray band observed with Fermi-hAT. Dur¬ 
ing a short time interval of the gamma-ray flare period 
( 0.2 days), there were simultaneous optical observations, 
whose results did not show any correlated variability 
with the gamma-ray flare as presented in H15. The X- 
ray observations in the period are available from Swift- 
BAT transient mon itor results Q by the Swift-BAT team 
(iKrimm et al.l 1201^ . The data provided an upper limit 
in the 15-50 keV band. We adopt the same values for Rq, 
F, T(jy, and U{jy as those in the previous section. By 
changing K', N^, and Bq, we attempt to fit the spectrum 
of the flaring period B in 2013. 

We consider one shell that contributes to the flare emis¬ 
sion. The energy diffusion coefficient and injection rate 
are slightly increased from the values in the steady model 
to K' = 1.3 X lO-'^ s-i (face = 1/(2A:') = 0.25IF7c), 
and = 2.5 x lO^® s”! « = 0.85{R/Rq)-^ cm”® s"!). 


^ http://swift.gsfc.nasa.gov/results/transients/weak/3C279/ 
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Fig. 1.— (a) Evolution of the electron energy distribution with increasing distance R in the steady model. The numbers beside each line 
denote R/Rq. The electron spectra during the acceleration phase {R < 2Rq) are plotted with blue lines, while those after the shutdown 
of the acceleration {R > 2Ro) are black, (b) P hoton spectrum for the active period in 2009. The open circles are measured flux points in 
period D (adapted from IHavashida et al.ll2012ll . The solid line is the model spectrum that is a superposition of the emissions from the all 
shells of /? > Rq. 


respectively. Hereafter, we call this the “fiducial” flare 
model. No significant concurrent flare in the optical 
bands implies that the optical photons are emitted from 
other steady components. As discussed in Section 4.3 in 
HI5, the lack of overall correlation between the optical 
and gamma-ray bands in 2013-2014 also suggests the dif¬ 
ferent origin of the optical component. The synchrotron 
flux of the flare should be below the observed flux level 
so that an upper limit for the magnetic field in the flare 
zone will be given. Here, we adopt a very low value of 
Ho = 0.25 G. 

Figure [2] (a) shows the evolution of the electron energy 
distribution in this flare model. Electrons are accelerated 
to higher energies compared to the case in the steady 
model. The power-law distributions above 10® eV are 
due to not only the larger K' but also the inefficiency of 
the synchrotron cooling owing to the low magnetic field. 
The secondary bumps at ~ 2 x 10® eV for R = 3-4i?o 
are attributed to the generation of secondary electron- 
positron pairs via internal 77 absorption. 

The resultant gamma-ray spectra shown in Figure [2] 
(b) agree well with the observed gamma-ray data. Here, 
we add an underlying component (that overlaps the solid 
gray line in the figure) to explain the gamma-ray flux be¬ 
fore the flare and the optical data. The flare spectrum 
has a higher synchrotron peak energy than the model 
in HI 5, since our flare model shows a drastic growth 
of the maximum energy of electrons compared to the 
steady model. As remarked above, the weak magnetic 
held strikingly increases the maximum energy of elec¬ 
trons. Even for this low magnetic held, the optical hux is 
slightly enhanced during the hare. The steady behavior 
of the optical light curve may prefer a weaker magnetic 
held, but we regard this as a conservative upper limit. 

The sharp cutoh at ^ 10^® eV in the photon spectrum 
is due to the 77 absorption inside the emission region. 
Some fraction of photons above 10^® eV escape from the 
shell. The model hux at the 100 GeV band is still higher 
than the detection limit for Cerenkov telescopes. How¬ 
ever, it should be noted that we have neglected the 77 
absorption after the escape from the shell. The absorp¬ 


tion by the broad emission lines during propagation may 
greatly suppress the hux around 100 GeV. 

Even for the same values of i?o and F as in the steady 
model, the light curve is well reproduced as shown in 
Figure [21 Thus, the emission zones of the intense hare 
in 2013 and the active period in 2009 may be located at 
similar distances from the central engine. The observed 
asymmetric prohle in the light curve is favorable for our 
simple one-shell emission-zone model. The strong cool¬ 
ing due to the external IC yields the rapid decay of the 
light curve. The evolutions of energy densities in Figure 
[4] clearly show the energy input by the SA and rapid cool¬ 
ing just after the end of the acceleration. At i? = 2i?o, 
the energy density ratio of the magnetic held to electrons 
is quite low as U'^/U'^ 4 x 10“®. 

The observational constraints, of course, do not deter¬ 
mine the model parameter uniquely. However, the essen¬ 
tial parameter for determining the gamma-ray spectral 
shape is only K' in our model (the role of N'^ is just nor¬ 
malizing the hux level, and the value of Bq does not affect 
the gamma-ray spectral shape). In Figure|21 we compare 
several photon spectral models derived with different pa¬ 
rameter sets. When we reduce K' by a factor of two 
from the hducial model (“Low- AT” model: K' if'x0.5, 
Ad —>■ Ag x4.8, the others are the same), the peak photon 
energy does not reach 10 GeV. Gonversely, we increase 
the diffusion coefficient as shown in the “High- if ” model 
(if' ^ if' X 1.5, Ni ^ N^x 0.038, Bq ^ Bq x 0.4, 
the others are the same). In this case, we need an even 
weaker magnetic held. The peak time of the light curve 
is delayed due to the lower compared to the hducial 
model. In Figure [21 we shift the light curve by 1.5 hr ear¬ 
lier. Other physical parameters (F, etc.) of the jet were 
derived from the steady model. However, as H15 sup¬ 
posed, we also try to increase F in our model. The initial 
radius should be increased as oc F^ to keep the variability 
timescale. Such an example (“High-F” model: F ^•rx2, 
Ro ^ RqX 4, ^ X 0.039, Bq ^ Bq x 0.16, if' 

is the same) is shown in Figure |21 Due to the relatively 
short face (= 0.181V'/c), the maximum electron energy 
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Fig. 2.— (a) Evolution of the electron energy distribution with increasing distance R in the fiducial flare model. The figure format is the 
same as that of Fig. 1(a). (b) Photon spectrum for the most intense flare in 2013. The open circles are measured flux points in period B 
(adapted from H15), and the 2cr upper limit for the hard X-ray is obtained from Swift-'QP^ data. The solid lines are the model spectra at 
observation times of 2 (black), 3.5 (blue), 6.2 (red), 11 (green), 20 (purple), and 35 (gray) hours. The model parameters for an underlying 
component, which is consistent with the gamma-ray flux the day before the flare event (0.19 x 10“^ photons cm~^ s“^ above 0.1 GeV) 
and the optical data, are the same as those of the steady model in ^ except for Bq = 3.8 G and A' = 7.3 x 10“^® s“^. The dashed lines 
show the flare components only. The observation time tobs is measured from the first arrival time of the photons escaped from the shell at 
R = Rq, Photons are supposed to be emitted from a single shell moving toward us. 


Photon Flux [Phs/cmVs] (>0.1GeV) 



Fig. 3.— Gamma-ray light curve of period B in 2013 with the 
flare model (solid black), and “High-F^” model (dashed orange). 
The original data were obtained from H15. In the model light 
curves, the underlying gamma-ray level is the same as in Fig. [2] 
(0.19 X 10“^ photons cm“^ s“^). 

grows as high as e' ~ 10^^ eV. A very low magnetic field 
{Ug/Ug < 6 X 10“®) is necessary again to suppress the 
synchrotron flux in the X-ray band. The strong cooling 
due to the higher oc makes the GeV spectrum 
too soft compared to the observed gamma-ray spectrum. 

5. DISCUSSION 

The simple SA model can reasonably explain the very 
hard spectrum and short variability in the intensive flare 
in 2013. The turbulence driving the particle acceleration 
may be generated by the hydrodynamical instability or 
the magnetic reconnection. 

Compared to the steady model for the active period in 
2009, the drastic alteration we need is the decrease of the 
magnetic field. The other parameters have almost sim¬ 
ilar values. The absence of the optical flare implies the 
weak magnetic field (< 0.25 G). The requirement of the 
magnetic field decrease at gamma-ray flare stages was 


Energy Density [erg/cm^] 



Fig. 4.— Evolution of the energy densities in the shell frame for 
the flare model. The solid black, blue, and red lines show the value 
of electrons, magnetic fields, and photons produced in the shell, 
respectively. The red dashed line denotes the value of the external 
photons. 

suggested by lAsano et al.l (|2014ll as well. The required 
low magnetic field seems irrelevant to the energy source 
for the particle acceleration. Therefore, a hydrodynami¬ 
cal instability is responsible for driving the SA. 

When r = 15, the variability timescale is consistent 
with Rq = 0.02 pc as shown in Figure [31 This distance 
from the engine also agrees with the constraint by the 
X-ray SSC component in the active period in 2009. The 
size of the central engine may be ~ 3rg ^ 8.6 x 10“® pc 
for the black hole mass of 3 x 10®Mq. If we adopt the sim¬ 
plest model for the jet acceleration d ue to the magnetic 
energy dissipation (iDrenkhahn 11200^ . the bulk Lorentz 
factor at i? = i?o should be < (i?o/3rg)®/® ^ 6, which is 
inconsistent with the postulated value of T. Given the 
variability timescale At, the initial radius can be scaled 
as Rq oc Atr^ when we change T. However, even in this 
case, the maximum Lorentz factor at Rq inferred from 
the magnetic dissipation model increases by a factor of 
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Fig. 5. — Comparison of the model spectra at fobs = 6.2 hr: the 
same model as in Fig. [2] (solid red), “High-i^” model (dashed or¬ 
ange), “Low-A"” model (dashed blue), and “High-F” model (dashed 
green). The steady underlying component is included as in Figure 
[2 The thin lines show the flare component only. 

only oc F^/^. For the jet acceleration model by the 
Poynting flux dissipation, not only the low magnetic field 
but also the short variability timescale are problematic. 
This problem is also raised for the very short gamma- 
ray flare (a few hundreds seconds) of B L Lac objects like 
PKS 2155-304 (|Aharonian et al.l 1200711 . For such SSC- 
dominant objects, however, the distance from the engine 
is not well constrained compared to FSRQs. 

The tiny change of the diffusion coefficient K', in spite 
of the drastic decrease of the magnetic field, seems enig¬ 
matic. The assumed value of g = 2 may be favorable 


for this invariant behavior of K'. In this low magnetic 
field case, the average energy gain per scattering may be 
proportional to /5 wj where /?w is the average turbulence 
velocity, rather than the Al fven velocity. The pitch an- 
gle diffusion approximation (iBlandford fc Eichler ifToS^ 
and power-law magnetic turbulence of 6B^ {k) = 6BQk~‘^, 
where k is the wavenumber, leads to AT oc 0^5BqB~^. 
The last factor of B~^ implies that electrons inter¬ 
act with higher (lower) amplitude turbulence at longer 
(shorter) wavelengths for a weaker (stronger) magnetic 
field. If 5Bq (X B, q = 2 results in AT oc /3|^, which 
is independent of B. Alternatively, magnetic bott les as 
“hard spheres” (jBeresnvak. Yan fc Lazarian Il20lil l may 
be formed in turbulence independently of the strength 
of the magnetic field. Such requirements for the turbu¬ 
lence property motivate us to probe the hydrodynamical 
instabilities in blazar jets. 
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